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ABSTRACT 

Algebraic  facts  about  the  space  of  polynomials  contained  in  the  span  of  integer  trans¬ 
lates  of  a  compactly  supported  function  are  derived  and  then  used  in  a  discussion  of  the 
various  quasi-interpolants  from  that  span. 
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SIGNIFICANCE  AND  EXPLANATION 


*  I 

\ 

The  linear  span  of  integer  translates  of  a  fixed  compactly  supported  function  pro¬ 
vides  a  particularly  simple  model  of  an  approximating  family  of  the  finite  element  type. 
The  approximating  power  of  such  a  span  (or,  more  precisely,  of  its  scaled  versions)  has 
been  known  for  some  time  to  be  characterizable  in  terms  of  the  space of  polynomials 
it  contains.  ■  '  (  •  ^  7 ' 

Recent  work  on  box  splines  has  provided  concrete  examples  of  interest  in  a  multivariate 
setting  and  so  rekin4ied  interest  in  the  space  The  report  derives  and  extends  specific 
information  about  ti^  contained  in  recent  work  by  Dahmen  &  Micchelli,  and  by  Chui, 
Diamond,  Jetter,  Lai  and  Ward,  but  does  so  without  reference  to  specific  properties  (such 
as  piecewise  polynomiality,  or  factorizability  of  the  Fourier  transform)  of  tp.  ^  •  . 

Understanding,  in  the  simplest  possible  and  most  efficient  terms,  of  the  approximation 
power  of  such  spaces  may  provide  the  necessary  insight  into  approximation  by  smooth 
piecewise  polynomials  on  regular,  and  perhaps  even  not  so  regular,  partitions. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary  lies 
with  MRC,  and  not  with  the  author  of  this  report. 
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THE  POLYNOMIALS  IN  THE  LINEAR  SPAN  OF  INTEGER  TRANSLATES 
OF  A  COMPACTLY  SUPPORTED  FUNCTION 

Carl  de  Boor 


This  note  was  stimulated  by  the  recent  papers  [CD85],  [CJW85],  and  |CL85]  in  which 
the  authors  take  a  new  look  at  the  space  of  integer  translates  of  box  splines  and,  in 
particular,  introduce  and  highlight  the  commutator  of  a  locally  supported  pp  function 
<p  of  several  variables.  The  intent  of  this  note  is  to  offer  alternative  proofs  of  some  of 
these  results,  and  to  point  to  some  connections  with  earlier  work  (e.g.,  [BH82/3],  [DM83], 
[BJ84]),  but  also  to  focus  more  attention  on  the  space  of  all  polynomials  contained  in 
the  span  of  the  integer  translates  of  the  box  spline  (or  other  compsurtly  supported) 

The  first  section  collects  simple  algebraic  facts  about  and  the  action  of  the  linear 

map 

-  f  ^  -»/(» 


The  second  section  records  that  is  invariant  under  differentiation  and  translation, 
and  brings  yet  another  characterization  of  jr^,  this  time  in  terms  of  the  Fourier  transform 
of  ip. 

The  final  section  makes  use  of  these  facts  about  in  a  discussion  of  the  various 
quasi-interpolants  available. 

Throughout,  I  will  use  standard  multi-index  notation.  I  find  it  convenient  to  use  the 
special  symbol  ||“  for  the  normalised  monomial  of  degree  a,  i.e.,  for  the  map  given 
by  the  rule 


With  this, 


:IR‘'  ^  IR  :  X  i“/a!. 


’Ta-  span(lf)^<^ 


denotes  the  space  of  all  polynomials  of  degree  <  a,  and 


TT*  :=  span(i|!^) 


\0\<k' 


7r<*  :=  span 


>!<*’ 


TT  :=  span! 


have  similarly  obvious  meaning. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-<X)41. 


1.  The  polynomials  Consider  the  span  of  integer  translates  of  a  compactly 
supported  function  ^  on  IR'*,  i.e., 

5  :=  :=  {(jp  *  c  :  c  e  (l.l) 

Here  1  use  the  convolution  product  notation 

IP  *  c  -  j)c{j)  (1.2) 


since  there  is  no  danger  of  confusion  with  either  the  continuous  or  the  discrete  convolution 
product.  I  find  it  convenient  to  use  the  special  notation 

<P*' f  •=<P*  ^  <p{- -  (1.3) 

iez*' 

in  case  /  is  a  function  on  51“^,  in  order  to  stress  the  semi-discrete  character  of  this  product. 
Further,  since  the  restriction  to  TL^  of  a  function  on  IR*^  occurs  often  here,  1  will  employ 
the  abbreviation 

/|  •=  /iz-' 

for  it. 

The  asymmetry  in  the  semi-discrete  convolution  product  (1.3)  is  not  all  that  strong 
since,  after  all, 

(p  **  f  z=  f  <p  on  TL^. 

This  implies,  e.g,,  that,  for  /  G  tt  (hence  f  ^  €  tt), 

(p  *' f  =  f  *' p  <=> 


hence 


TT  ^  •  — 

{f  e  n  :  p  *'  f  e  n}  =  {f  e  n  :  p  *'  f  =  f  *'  p}. 

(1.4) 

It  also  implies  that 

p  *'  f  =  f  *'  p  for  all  f  E  S, 

(1.5) 

since,  for  /  =  *  c. 

p  *'  f  =  p  *  {pf  *  c) 

=  <p  *  (c  *  v?|) 

—  {p  *  c)  *  p\  =  f  *'  p. 

As  a  consequence,  one  gets  the  inclusion 

TT  n  s  c  {f  e  TT  :  p  f  -  f  *'  p} 

(1.6) 

and  the  conclusion  that 

p,'-.f^p,'f 

2 


maps  into  tt  P  5.  This  implies  that  there  must  be  equality  throughout  (1.6)  as  soon  as 
the  linear  map 

can  be  shown  to  be  1-1 .  But  that  is  easy  to  do  under  the  assumption  that  v?  is  normalized, 


51  =  1- 


For,  under  this  zissumption, 


for/€7r^,  <p*'/  =  “  EjeZ-*  ~ 


^<deg/ 


since,  for  each  j,  /-/(•-  j)  G  7r<aeg/- 


The  salient  facts  of  this  discussion  are  gathered  in  the  following. 
Proposition  1.  //v^  is  normalized,  then 

^<p‘  =  {f  C  :  <p  *' f  £  n)  =  {f  e  TT  :  ip  *' f  =  f  *' <p} 


7r  n  5 


Further,  L  ;=  is  onto,  and 


=  {/ e  TT  :  V?  *' /  6  /  +  7r<deg/}. 


U  ;=  1  -  L 


is  degree-reducing.  In  particular, 


Ip  O  TTq)  — 


(1.10) 


As  a  consequence,  =  0  on 


TTyj .  TT^  P  ^<fc- 


Therefore 


(v.)  '  =  - 


(1.11) 


Note  that  is  necessarily  finite  dimensional,  since  p  is  compactly  supported. 
Precisely,  for  any  bounded  set  D.  the  set 

A(D)  {q  t  :  p(-  -  a)  0} 


is  finite,  hence  if  D  also  has  interior,  then 


dim  TT-r  ^  dim  ^  AID)  '  oc. 


The  sharpest  bound  attainable  this  way  for  a  piecewise  continuous  ip  would  be 


dim;r^  <  max  #y4({x}).  (l-U) 

Z 

In  any  case,  this  implies  that 

X'*  =  1  +  +  (/*  +  ••■  , 

with  the  Neumann  series  actually  finite. 

The  assumption  that  p  be  normalized  is  no  real  restriction  except  when 

V3(i)  =  0. 

]€Z‘‘ 

In  this  case,  (1.7)  shows  L  to  be  degree-reducing,  hence  in  particular,  not  invertible. 
Consequently,  tt  n  5  may  be  strictly  smaller  than  For  example,  with  p  =  1  on  [-1,0[, 
=  -1  on  [0,  l{,  and  =  0  otherwise,  =  ttq  ^  {0}  =  tt  n  S. 


r»,  w'.  ir.TTw  ■ 


2.  Invariance  Denote  by  E  the  multivariate  shift,  i.e.. 


£“/  ;=  /(•  ^  a),  a  e  5Z  . 

While  it  is  obvious  that  commutes  with  E,  hence  TTip  is  invariant  under  E,  some  of  the 
other  properties  of  derivable  from  this  fact  may  not  be  as  immediate. 

Proposition  2.1  The  linear  map  L  =  commutes  with  differentiation,  hence 

with  translation,  i.e., 

=  D'"!,  Vae  2Z+^  E^L  ^  LE\  e  IR‘‘.  (2.1) 


Proof  Since  is  a  finite-dimensional  polynomial  subspace,  there  exists,  for  each 
a  6  a  weight  sequence  w  of  finite  support  so  that 

=  '^  w(/3)E^  on  TT^.  (2.2) 

(E.g.,  with  t,  the  Lagrange  polynomials  for  the  points  0, . . .  ,fc  :=  max  deg  tt^,  we  have 

p  =  Y.  e'’£^p(o) 

O<0U)<k 

for  all  p  e  trjt(IR)  0  •  •  •  (gi  7r;t(IR)  D  hence  w{0)  :=  D“£^(0),  all  /?,  would  do.)  Thus, 
LE  -  EL  implies  LD  =  DL.  But  this  finishes  the  proof  since 


Remark  The  argument  shows  that  any  ^-invariant  polynomial  subspace  is  D- 
invarianl,  hence  even  translation-invariant,  i.e.,  for  any  linear  subspace  P  of  tt. 


Vq  e  E^P  C  P 


Va  €  TL‘^  D'^P  C  P 
Vy  e  E^P  c  p 


('orollary  is  D-invariant  and  translation-invariant. 


As  a  simple  consequence,  consider  the  polynomials  g^,  defined  in  jCJW85j  by  the 


recurrence 


-  Y  ‘’s-’ 


and  then  shown  to  satisfy 


(2.6) 


=  Y1  9a{j)^{x-j) 

in  case  'al  <  m  and  7r<m  C  In  other  words,  Therefore,  on  first  reading, 

I  thought  that  the  recurrence  relation  (2.5)  was  a  consequence  of  the  fact  that  L  is  “unit 
upper  triangular”.  In  fact,  the  recurrence  can  be  derived  from  the  identity  DL  =  LD. 

For  this,  recall  that  the  Appell  sequence  for  a  continuous  linear  functional  n  on 
C(IR  )  with  p(l)  =  1  is,  by  definition,  the  sequence  (g^)  determined  by  the  conditions 

Qa  ^  ^ori 

There  is,  in  fact,  exactly  one  such  sequence  for  given  /i  since  the  linear  system 

„£)'>(  r a,)  = 

7<a 

for  the  power  coefficients  (a.,)  for  has  a  unit  triangular  coefficient  matrix.  Backsusbsti- 
tution  therefore  provides  the  formula 

=  ir  -  X]  (2-5') 

whose  correctness  can  also  be  verified  directly  by  induction  on  a: 

=  Mlf-"  -  Mir-"  =  0 

for <  a,  while  =  m(1)  =  1-  With  existence  and  uniqueness  established, 

facts  about  the  Appell  sequence,  such  as  symmetries  which  reflect  those  of  n,  or  that 
D^ga  =  9a -0^  follow  immediately. 

In  our  case,  n  :  f  tp  *'  f(0],  hence,  for  .  tt  ,. 

=  fiD^g^  ■=  p  *'  (£>‘V)(0)  iy\p  *'9a){0), 

which,  together  with  the  fact  that  p  *'  ga  6  L-n^  -  shows  that 

(2.6') 

The  resulting  different  normalization  of  as  compared  with  (2.5)  avoids  all  those  facto¬ 
rials. 


Dahmen  and  Micchelli  DM83  consider  the  polynomial  space 

{p  t  TT  :  p[D)p  -  0  on  27rffi'^\0}. 


(2.7) 


with  the  Fourier  transform  of  if.  It  seems  slightly  more  convenient  to  consider  instead 

;=  {p  G  7r  ;  p(-iD)ip  =  0  on  27rZZ‘^\0}. 


They  prove  that  any  affinely  invariant  (i.e.,  translation-  and  scale-invariant)  subspace  of 
(2.7),  hence  of  n<p,  is  contained  in  n^.  But  their  proof  can  be  made  to  show  more. 

Proposition  2.2  js  the  largest  E-invariant  subspace  of  Y\^. 

Proof  The  proof  in  [DMSS]  is  based  on  the  observation  that,  by  Poisson’s  sum¬ 
mation  formula, 


f  *'  p(i)  =  “  «)p{«)  =  ^0(27rQ), 

a  a  a 

while,  for  any  p  €  tt,  the  function  tp  :  y  <p(x  -  y)p(y)  has  the  Fourier  transform 

ip(y}  =  e~'^^(^p(x  -  iDjf'ji-y). 

If  now  p  &  P,  with  P  an  ^-invariant  (hence  D-invariant)  subspace  of  11^,  then 
ip(2na)  =  (^p(x  -  iD)ip^{2na)  =  p[~iD)ip^[2‘no)  =  0 


for  a  0,  hence 


■p  =  (p(l  -  iDj^)(0) 

=  5;oWS-.i?!V(0) 

a 

=  p{x]f[0)  -  D'^p{x)[  -  ?£>|“v^(0). 

\a'>0 


showing  that  f  *'  p  iz  tt,  i.e.,  p  €  tt^. 

On  the  other  hand,  if  p  £  then 

f*'p=  -  iD)yf(-2na) 

a 

is  a  polynomial,  and  this  is  possible  only  if 

p(-  -  iD)f{2na)  =  0  va  0. 

showing  that  p  £  fl^. 

Corollary  tt,  ^  :  rr^  arid-.  ^  -  k^h  ' 


Proof  If  i  a  +  01  <  k  +  h,  then  either  |a|  <  k  or  else  10\  <  h,  hence  i'll  <  k  +  h 
implies 

1  -  *  V')(27ry)  =  ^  l-iDf<p{2Tij)l-iDfrl>{2nj)  =  0 

for  j  6  E^O. 

While  TT^  has  been  shown  in  1BH82/31  to  be  scale-invariant  in  case  v?  is  a  box  spline, 
it  is  not  clear  that  tt^  is  necessarily  scale-invariant  for  arbitrary  (p.  For  this,  I  note  that  a 
polynomial  subspace  P  is  scale-invariant  if  and  only  if  P  stratifies,  i.e.,  P  =  P 
with 

span 

Hence,  span{j|^’°  -i-  1}  provides  a  simple  example  of  an  F^-invariant  polynomial 

subspace  which  is  not  scale-invariant. 


3.  Quasi-interpolants  The  space  is  of  interest  because  it  characterizes  the 
local  approximation  order  obtainable  from  5,  or,  more  precisely,  from  the  scale  (5/,) 
associated  with  S.  To  recall, 

Sh  ■=  Oh(S), 

with 

Ohf  f{x/h.). 

Further,  the  local  approximation  order  of  S  is  the  largest  k  for  which 

dist{f,Sh)  =  0{h’^) 

for  all  smooth  /,  with  the  distance  measured  in  some  norm,  e.g.,  the  max-norm  on  some 
bounded  domain,,  and  the  support  of  the  approximation  to  /  within  h  of  the  support  of 
/• 

In  iFS69j,  Fix  and  Strang  give  a  characterization  of  the  local  approximation  order 
from  the  scale  {5^.)  which,  in  the  terms  of  Section  1,  can  be  phrased  thus:  it  is  the  largest 
k  for  which 

U  :=  1  -  <p  *'  is  degree-reducing  on  7r<fc.  (3.1) 

Proposition  1  shows  that  we  can  state  this  condition  more  simply  as 

(3.2) 

To  be  precise,  'FS69]  consider  the  “controlled”  approximation  order,  which  turns  out  to 
be  the  same  as  the  local  approximation  order;  cf.  [BJ84]. 

Fix  and  Strang  use  in  their  proof  a  quasi-interpolant  whose  construction  relies  on 
Fourier  transform  arguments  which,  in  a  univariate  context,  can  already  be  found  in 
Schoenberg's  basic  spline  paper  [846]  and  which  appear  in  the  proof  of  Proposition  2.2. 
This  makes  it  easy  to  recall  their  construction  here. 

Define  the  quasi-interpolant  Q  on  tt  by  the  rule 

Qf 

with 

Ff  :=^a^(-iDrf 

a 

and  ,3)(0)  the  Taylor  coefficients  for  1/ip.  Dahmen  and  Micchelli  [DMSS 

prove  that  Q  reproduces  any  affinely  invariant  subspace  of  (2.7),  but,  again,  their  argument 
supports  a  stronger  claim,  viz.  that 


s' 

V 


For,  if  p  6  Ttip,  then  also  Fp  6  since  is  D-invariant;  hence,  by  (2.8), 

Qp=5^(Z?^Fp)i-.£>iV(0) 

Or 

=  i£>r^(o) 

Q  ^ 

^Y.i-'oyp  E  iDf(i/,5)(o)|Dr,p(o) 

1  a+0=-i 

=  Y^i-iDyp  6o^  =  p. 

7 

The  construction  is  finished  by  noting  that  (3.3)  only  depends  on  the  action  of  F  on  tt^, 
hence  a  local  quasi-interpolant  on  smooth  functions  which  reproduces  tt^  can  be  obtained 
in  the  form 

Q/;=7r^*'(A*/),  (3.4) 

with 

(A*/)(x):=  A/(.  +  i),  (3.5) 

and  A  any  locally  supported  linear  functional  which  agrees  on  with  p  *-+  Fp(0). 

The  construction  idea  in  [BH]  seems  more  direct:  There  the  locally  supported  bounded 
linear  functional  (on  whatever  normed  linear  space  X  you  may  wish  to  carry  out  approxi¬ 
mation  from  S  n  Af)  is  constructed  as  an  extension  of  the  linear  functional 

p^  {L~^p){0).  (3.6) 

Since  L  =  commutes  with  E,  so  does  L'*.  Thus,  for  p  G  tt^. 


hence 


(L  ‘p)(i)  ^  (L  ^p(--t- j)){0)  =  (A*p)(i), 


Qp-  (L  *p)  =  P- 


In  order  to  obtain  a  quasi-interpolant  of  the  optimal  order  k,  the  extension  A  only 
needs  to  match  (3.6)  on  7r<fc.  For  example,  one  obtains  the  Strang-Fix  quasi-interpolant 
by  expressing  the  extension  as  a  linear  combination  of  the  linear  functionals 

f  ^  {-iD)^f(0),  a  <k,  (3.7) 


A/=  E  aa(  «X>)V(0). 


i.e.,  in  the  form 


The  weights  Oq  are  uniquely  determined  by  the  requirement  that  this  linear  functional 
match  (3.6)  since  (3.7)  is  maximally  linearly  independent  over  7r<jt.  In  particular, 

a,  =  L-‘It-r(0)  =  i“(,,(0), 

by  (2.6’).  This  shows,  incidentally,  that 

lDr(l/^)(0)  =  t-g<,(0). 

If  point  evaluation  is  continuous  on  X,  then  the  linear  functional  A  can  be  written  zis 
a  linear  combination  of  evaluations  at  integer  points  near  0.  For,  by  (1.11), 

=  (1  +  t/  +  •••  + 

while,  from  (1.9), 

(y/)0)  =  (c '/;)(». 

with 

c  -  ^1 

and  i  the  unit  sequence,  i.e.,  6{j)  =  6j.  Hence 

(t-‘p)(0)  =  pl*'(0) 

with  obtained  inductively  in  the  following  computation: 

f  0,  if  r  =  0; 

pi'!  :=  (3.8) 

P|  +  c  *  pl’'  *1,  if  r  >  0. 

This  gives 

(Z/~V)(0)  ==  P&7r<k 

jez' 

with  the  weight  sequence  C  of  finite  support  since  c  has  finite  support. 

This  construction  was  arrived  at  by  different  means  by  Chui  and  Diamond  [CD85], 
who  added  the  following  very  useful  observation.  If  v?  is  symmetric,  then  U  reduces 
the  degree  by  at  least  2,  since  (1.7)  can  then  be  written  in  the  form 

for  /  €  7r^,  Ip  *'  /  =  /  51  ^  -  2/  +  /(•  -  j))<p(j)/2.  (1.7') 

This  implies  that,  on  already  t'  *  vanishes,  hence  only  half  the  iteration  (3.8)  is 
necessary  in  this  case. 

Even  for  a  symmetric  ip.  the  support  of  the  resulting  A  may  be  far  from  minimal. 
Since  we  are  only  interested  in  extending  a  linear  functional  from  a  support  consisting 


of  (dim  TT^p)  points  is  sufficient.  These  points  can  be  chosen  from  Tu  since  TL^  is  total  for 
TT.  It  would  be  interesting  to  find  out  whether  they  could  be  chosen  as  neighbors. 

Such  questions  of  minimal  support  for  A  have  been  answered  quite  elegantly  by  Dah- 
men  and  Micchelli  in  case  ^  is  a  box  spline.  They  find  in  [OMSS]  that  the  (dim 
integer  points  in  the  (right-continuous)  support  of  v?  are  linearly  independent  over  )r^, 
and  so  conclude  the  existence  of  an  extension  from  involving  just  these  (dim  n^)  point 
evaluations. 

1  note  that  the  quasi-interpolant  construction  in  {BJ84|  takes  the  opposite  tack.  In¬ 
stead  of  constructing  an  appropriate  A  as  a  linear  combination  of  certain  point  evaluations, 
a  compactly  supported  function  e  S  is  constructed  there  so  that  already  reproduces 
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